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Abstract 

A new technique for measuring dynamic gas 
velocity and temperature is described. The technique is 
based on molecular Rayleigh scattering of laser light, so 
no seeding of the flow is necessary. The Rayleigh 
scattered light is filtered with a fixed cavity, planar 
mirror Fabry-Perot interferometer. A minimum number 
of photodetectors were used in order to allow the high 
data acquisition rate needed for dynamic measurements. 
One photomultiplier tube (PMT) was used to measure 
the total Rayleigh scattering, which is proportional to 
the gas density. Two additional PMTs were used to 
detect light that passes through two apertures in a mask 
located in the interferometer fringe plane. An 
uncertainty analysis was used to select the optimum 
aperture parameters and to predict the measurement 
uncertainty due to photon shot-noise. Results of an 
experiment to measure the velocity of a subsonic free 
jet are presented. 

Introduction 

A variety of nonintrusive, laser based flow 
diagnostics (Rayleigh scattering, LDV, PIV, EXJV, LIF) 
are routinely being applied for time average and 
instantaneous planar measurements of velocity, density, 
temperature, and species concentrations. However, an 
important class of measurements is not addressed by 
current laser diagnostics. High frequency response 
dynamic measurements of flow parameters cannot be 
made with existing techniques. This type of time 
history data is needed to determine, for example, 
density and velocity spectra, density-velocity 
correlations, and two-point correlations. Although 
Laser Doppler Velocimetry (LDV) can typically 
achieve data rates of a few tens of kHz, it is difficult to 
achieve much higher rates without introducing larger 
amounts of seed material into the flow. In practice, 
LDV normally provides mean velocity and a measure of 
turbulence intensity. Furthermore, the measurements 
are random in time, which makes it difficult to obtain 
time history data needed for power spectra. In turbulent 


flows, LDV measurements are beset by a variety of so 
called biasing errors, caused by correlations between 
the measurement rate and flow properties. Planar 
techniques, such as PIV, provide a large number of 
simultaneous measurements in the measurement plane, 
but are generally limited to low sampling rates, 
determined by the pulse repetition rate of the laser and 
by the time needed to transfer image data from the 
camera. 

New techniques are needed to provide 
nonintrusive, dynamic measurements that can provide 
data similar to that provided by hot wire anemometers. 
Since it unrealistic to expect (at least at the present 
time) to make measurements at a large number of 
locations and at a high sampling rate, it seems prudent 
to develop a laser diagnostic capable of point 
measurements at a high sampling rate. Such time 
resolved point measurements are in a form familiar to 
fluid dynamics researchers. 

In this paper, we describe a laser diagnostic for 
dynamic measurement of gas properties based on 
molecular Rayleigh scattering. Because the measurement 
is based on molecular scattering, rather than particle 
scattering, no seed material need be injected into the flow 
under study. One difficulty with LDV and PIV is that 
they require the flow to be "seeded" with micron sized 
seed particles to provide a sufficient concentration of 
scattering centers. The reliance on seed particles presents 
a number of difficulties. It is difficult to inject a uniform 
cloud of seed into many test facilities. And the seed 
material must withstand the flow environment. For high 
temperatures, refractory seed materials are necessary. In 
addition, seed material can contaminate the facility, 
coating surfaces and windows. Facility engineers are 
often wary of the introduction of this foreign, often 
abrasive material, into their equipment. A further 
limitation of particle scattering methods is that the 
particles, although small, may not be able to follow large 
flow accelerations, thus introducing inaccuracy in the 
prediction of the gas velocity. 
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These problems associated with particle 
scattering measurements are eliminated if molecular 
scattering is used, since the gas molecules that constitute 
the flow under study are used as the scattering centers. 
The simplest molecular scattering based diagnostic is 
Rayleigh scattering. The frequency spectrum of Rayleigh 
scattering is closely related to the velocity distribution of 
the scattering gas. The spectrum may be analyzed to 
determine temperature, density, and velocity. Density is 
simply proportional to the total scattered light; 
temperature is obtained from the width of the Rayleigh 
spectrum; and one component of velocity is proportional 
to the shift of the Rayleigh spectral peak from the 
frequency of the incident light. Because the spectral 
width is also a function of the molecular weight of the 
gases in the flow, knowledge of the gas composition is 
generally required. This, however, is not a concern in the 
proposed work, which is directed toward air flows where 
the composition is well defined. However, because 
velocity is determined from the peak of spectrum, it is 
independent of the gas composition. Rayleigh scattering 
is particularly suitable for measurement of supersonic and 
hypersonic velocities, where the mean molecular velocity 
(flow velocity) is greater than the random molecular 
velocity (temperature). One quantity of particular 
interest in the study of aerodynamically generated sound 
that can be obtained from Rayleigh scattering 
measurements is pV (the product of gas density and 
velocity). 

The key to obtaining high frequency response 
data is to combine a high power, continuous wave laser, 
or a high repetition rate pulsed laser, with high efficiency 
collection optics. Density measurements are relatively 
straightforward, since the Rayleigh scattered light is 
directly proportional to the density. Velocity and 
temperature measurements are more challenging because 
they must be determined from the spectral characteristics 
of the Rayleigh scattered light. 

Both atomic/molecular vapor absorption filters 
and Fabry-Perot interferometers can be used to measure 
spectral characteristics that are related to the gas 
velocity and temperature. Iodine vapor is commonly 
used as an absorption medium because it has a rich 
absorption spectrum in the green wavelength region ^ 
Use of absorption filters requires a tunable laser that 
can be tuned to a suitable absorption line. One 
approach for obtaining real-time velocity measurements 
is to use an absorption filter to analyze frequency 
modulated Rayleigh scattered light" A Faraday cell 
has been used as a frequency discriminator for seeded 
flows'^. Fabry-Perot interferometers, on the other hand, 
act as bandpass filters. They are quite flexible in that 
they can be used at any laser frequency, and the width 


of their passband can be easily adjusted. However, the 
transmission is also a function of the angle of the light 
beam in the Fabry-Perot interferometer cavity relative 
to the optical axis. Fabry-Perot etalons have been 
proposed for analyzing Rayleigh scattering from the 
atmosphere to obtain wind velocity and temperalure\ 

The approach taken in the present work is to 
use a Fabry-Perot interferometer with fixed mirror 
spacing. Characteristics of the transmitted light are 
measured using a minimum number of fast photo 
detectors. This allows a high sampling rate to achieve 
the desired dynamic response. 

We first establish the lower bounds for 
uncertainty in the measurements of velocity and 
temperature caused by photon statistical noise. We then 
examine a particular scheme based on a Fabry-Perot 
interferometer and three photomultipler tubes and 
determine the lower bounds for measurement errors. A 
similar approach has previously been used to determine 
lower bounds for temperature and velocity uncertainties 
for measurements obtained with a scanning Fabry-Perot 
interferometer^, as well as for an imaging system^. The 
results of a preliminary experiment to evaluate this 
technique are then presented. 


Theory 

For the purposes of this paper we consider 
Rayleigh scattering from a low density gas, which has a 
Gaussian spectrum given by 


S(f-fjdf=-^exp- 
aK 


aK 


df 


( 1 ) 


where /o is the laser frequency and u is the mean gas 
velocity. Note that the spectral peak is shifted by a 
frequency proportional to the component of the bulk 
velocity in the K direction. The spectral width is 
proportional to the square root of the gas temperature. 
It is convenient to introduce the velocity component 
= K* utK^ which represents the measured velocity 
component. The interaction wave vector is K = kj-ko 
(with ko and ks being the wave vectors of the incident 
and scattered light), and a = (2x77m)'^‘ is the most 
probable molecular speed (with k being Boltzmann's 
constant, m the molecular mass, and T the gas 
temperature). (Note that a = (2/)^^^ q, where is the 
speed of sound and yis the ratio of specific heats.). 

The assumption of a Gaussian shaped Rayleigh 
scattering spectrum is only valid if 


where p is the gas pressure and X] is the shear viscosity. 
For higher density gases (y ~ 1), a more detailed kinetic 
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theory model, such as the Tenti S6 model^ is required to 
describe the Rayleigh scattering spectrum. To simplify 
the analysis presented here, we use the Gaussian 
spectral model. However, as shown in Figure 1, the 
Tenti spectrum differs significantly from the Gaussian 
spectrum for the conditions encountered in this work ( v 
= 0.79). A more accurate analysis would use the Tenti 
spectral model. 

Since Rayleigh scattering is a relatively weak 
process, the uncertainty in the measurements often is set 
by the photon statistical noise (shot noise), which 
determines the lower bound on measurement 
uncertainty. For example, the variance in the number of 
photoelectron counts for a Poisson process is equal to 
the mean number of counts. Thus the lower bound for 
the relative uncertainty in the measurement of gas 
density p, is equal to the square root of the variance 
divided by the mean number of counts. The relative 
uncertainties can be written^^ 


(p) 

P 



crK)= 


(2) 


(2N«) 


These relations provide a lower bound for the 
measurement uncertainties, which can only be 
approached if the shot noise is the dominant noise and if 
an ideal instrument is used. Note that the lower bounds 
for the uncertainties of density, temperature, and 
velocity are all inversely proportional to the square root 
of the number of detected photons. 

Estimates of the measurement uncertainty for 
the technique described here, where the Rayleigh 
scattered light is analyzed with a planar mirror Fabry- 
Perot interferometer, are obtained by numerically 
calculating the Cramer-Rao lower bound*®. The 
variance of the estimate of a parameter a, (e.g., 
temperature or velocity) is given by 


V(a,) = [r'l, 


where, for Poisson statistics, T is the Fisher information 
matrix with elements 


n-I 




a Nn)d{N, 


da- 


da. 


(4) 


and the summation is over the number of 
photodetectors. For example, in previous work** using 
a CCD array detector, the summation would be over the 
pixels in the image being analyzed. For this work, 
however, a CCD camera cannot be used because of its 
slow readout rale (typically on the order of 1 second for 


scientific grade cameras). To achieve the high data 
acquisition rate desired for dynamic measurements, we 
use photomultipliers tubes (PMTs), which combine 
reasonable quantum efficiency with high gain and low- 
noise. One PMT is used to measure the total Rayleigh 
scattered light (i.e., the spectrally integrated light), 
which gives a measurement of the gas density. Two 
additional PMT’s are used to detect the Rayleigh 
scattered light that is transmitted through a planar 
mirror interferometer. 

Consider an experiment with the following 
parameters: a laser (532 nm) with output power = 4 
W, air at STP, probe volume length L, = 1 mm, f/4 
collecting optics (i.e., solid collection angle Q = 0.05 
sr), and efficiency factor £= 5 The rate of detected 
photons, given by 


N,= 


€ P nL X il( dcT ^ ^ 

sm~ X 


he 


da; 


= 39 million counts / sec 


(5) 


In this equation, n is the gas number density, da/da is 
the differential scattering cross section, x is the angle 
between the electric field vector of the (linearly 
polarized) incident light and the direction of the 
scattered light, h is Planck’s constant, and c is the 
velocity of light. If we wish to sample at a 10 kHz rate, 
the total number of detectable photons in each sampling 
period would be 3900. Equation 2 gives the lower 
bound for measurement of density, temperature and 
velocity for each sample period: 

= L6 %, = 23% , a(uf : ) = 3.6 m/sec 

P T 

This shows that high sampling rate Rayleigh scattering 
measurements are at least feasible. It must be 

emphasized, however, that these values represent the 
best possible measurements given an ideal instrument. 
In practice, we are limited to instruments, such as the 
Fabry-Perot interferometer used here, which result in 
significantly higher uncertainties, as described below. 
Achievement of the uncertainties given above will 
require a photoelectron detection rate considerably 
above those given in this example. This can be 

achieved by using higher laser power, larger collection 
f-number, or higher quantum efficiency photodetectors. 

Optical configuration 

In order to perform an uncertainty analysis, we 
first describe the setup (fig. 2) used for the experimental 
work. Light from a 5W, 532 nm, single-frequency. 
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Nd: Vanadate CW laser is focused by lens LI (300 mm 
focal length) to a 150 |itm diameter beam at the probe 
volume. The 90^ Rayleigh scattered light is collected 
and collimated by lens L2 (250 mm focal length) and is 
then split into two paths with an uncoated optical flat. 
About 10% of the light is reflected and focused by lens 
L3 onto PMT 1 (type R3896). This detected signal is 
proportional to gas density. Light in the second path is 
directed through a planar mirror Fabry-Perot 
interferometer (70 mm dia. mirrors, 90 % reflectivity, 
15 GHz free spectral range). The light exiting the 
interferometer is focused by the fringe forming lens, L4 
(135 mm focal length, f/2). An addition lens L5 (50 
mm focal length, f/1.2) is used to form a magnified 
image of the fringe pattern. The effective focal length 
of combination of lenses L4 and L5 is 1823 mm. The 
fringe pattern formed at the focal plane of this lens pair 
is analyzed to determine the velocity and temperature. 
Typical calculated images of the inner fringe of 
Rayleigh scattered light from a gas with zero bulk 
velocity and a temperature of 3(X) K are shown in figure 
3. One image shows Rayleigh scattered light from an 
extended region, which results in the full circular fringe; 
the other image shows the image resulting from a line 
source as was used in the experiment. 

Because, in the experiment, the laser beam is 
focused to a line in the probe volume, the fringe exists 
only along that line. Also, to simplify the setup and to 
use a minimum number of photodetectors, only one half 
of the fringe pattern is utilized. A mask having two 
apertures, shown as white rectangles in figure 3 and 
referred to as the “inner” and “outer” regions, is located 
in the fringe plane. Light passing through the inner and 
outer regions is directed to two PMT’s (type 8645). As 
the velocity and temperature of the gas change, the 
relative distribution of the light detected by the two 
PMT's changes. Note that a flow in the direction 
opposite the K vector results in a negative frequency 
shift. Use of additional detectors would improve the 
measurement uncertainty, but at an increase in the 
system complexity. 

Additional optics were included to provide for 
a reference image consisting of light at the unshifted 
laser frequency. To accomplish this, several 
components could be placed in the optical path. When 
placed in the beam, mirror M directed laser light 
through neutral density filter NDF onto a diffuser 
assembly that could be inserted at the measurement 
location. This assembly consisted of a diffuser that 
scattered light into the collection optics and an iris to 
limit the spatial extent of the light. Finally, a prism 
assembly could be placed in the light path between the 
Fabry-Perot interferometer and the PMT’s. This served 


to direct light into a standard video camera. The video 
signal from this camera was digitized by a frame 
grabber card in a 486 PC. A computer program 
analyzed this image and generated signals to control the 
Fabry-Perot mirror alignment. The mirror M, diffuser 
assembly, and prism assembly were mounted on 
remotely controlled pneumatic linear actuators. 

Uncertainty analysis 

We now calculate the Cramer-Rao lower 
bounds for velocity and temperature measurements 
based on the experimental setup described above. This 
allows us to conduct parametric studies to determine the 
optimum configuration for Fabr)^-Perot interferometer 
and for the light detection system. The observed fringe 
pattern for Rayleigh scattered light can be expressed as 

<N^> = jjA,S,{f)l,,{f,e)dfdQ ( 6 ) 


where <Np^> is the expected number of detected 
photons in the area. A/? is the number of photons that 
would be detected if the Fabry-Perot interferometer 
were removed from the system, and Ipp is the Fabry- 
Perot instrument function given by^“ 


I fP {y/) - 


1 + Fsin^ 



(7) 


where y/ is the phase change (neglecting any phase 
change on reflection) of the light between successive 
reflections given by 




f jji dcosO^ 

c 


( 8 ) 


Here, jj is the refractive index of the medium in the 
Fabry-Perot cavity, d is the Fabry-Perot mirror spacing, 
$r is the angle between the ray and the optic axis, and F 
= l/(sin‘(7i/2/V£) where Np is the effective finesse. In 
general, the image of a monochromatic extended source 
located in the object plane consists of a series of 
unequally spaced concentric rings. In this work, 
however, the field of view is restricted to include only 
the inner fringe as shown in figure 3. 

We now calculate the lower bounds for 
velocity and temperature measurement uncertainties 
based on the detected light from the inner and outer 
regions. The uncertainties are evaluated as a function 
of two parameters and np which characterize the two 
active regions in the fringe plane. Because the number 
of photoelectron counts from the two regions is a 
function of both temperature and velocity, a 
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measurement of the photoelectron counts from a single 
region cannot be used to determine velocity. 

For a given wavelength the order of interference at the 
center of the fringe pattern is given by = (/? + n^)l A,> , 
where n and represent the integral and fractional parts 
of the order, respectively^ I A bright fringe is located on 
the optical axis if n, = 0. Off axis, the order of 
interference at angle 6 is given by 2^ cosd = (n + — 

/ 7 ^)/Ao, where fu represents the change in the order of 
interference between optical axis and the angular radius 
8 . (Note that a bright fringe occurs at )• For 

small angles, = 1^0 '/A^ . The spectral interval 
included in the region from the optical axis out to 6 is 

where the free spectral range FSR = cP.d, It 
follows that the spectral interval between orders and 
nia is {Uk 2 ~ )FSR, For the purpose of analysis, it is 
convenient to use the fractional order to describe 
locations in the Fabry-Perot fringe plane. Because we are 
working in air, we assume that /i =1 in the remainder of 
this analysis. 

We can now define the boundaries of the 
fringe plane mask, as shown in figure 3, in terms of 
fractional orders: Ho is the fractional order of the fringe 
corresponding to the laser wavelength; (no - n/) and (no 
+ Hf) correspond to the inner boundary and outer 
boundaries of the central stop, respectively; and n^ 
corresponds to the outer boundary of the aperture. 
(Here we assume that the central stop of the mask is 
located at fringe position.) Note that the actual radius 
of a fractional order n^ in the focal plane of the fringe 
forming lens is given by the product the angle = (^o 
njd and the focal length of fringe forming lens. 
Also note that the image of the fringe pattern projected 
back to the probe volume is similarly given by the 
product of the angle and focal length of the collimating 
lens. For the mask used in this work, the active area of 
the probe volume is about 1 mm long. 

The lower bounds for velocity and temperature 
were numerically evaluated as a function of the 
parameters and «f, which characterize the dimensions 
of the two active areas. The results are shown in figure 
4. These plots are for a particular velocity (0 m/sec) 
and temperature (300 K). The order of the outer 
boundary of the outer aperture was = 0.341. Note 
that the minimum uncertainties for velocity and 
temperature occur at different values of n^ and «f. It can 
be seen that the values of no and that result in a 
minimum velocity uncertainty (about 13 m/sec) are 
0.087 and 0.027, respectively. On the other hand, the 
minimum temperature uncertainty (about 34 K) 
corresponds to «o and n^ values of 0.14 and 0.08, 
respectively. Other temperatures and velocities would 
result in other optimum values of «o and «f. For this 


study, we chose the configuration that results in the 
minimum velocity uncertainty at the expense of higher 
temperature uncertainty. We show in figure 5 the 
velocity and temperature uncertainties as a function of 
velocity and temperature for no = 0.087 and /if = 0,039 
(close to values that result in minimum velocity 
uncertainty at 0 m/sec and 300 K). The minimum 
velocity uncertainty is about a factor of four greater 
than that calculated above for an ideal instrument, and 
the minimum temperature uncertainty is about a factor 
of five greater than that of an ideal instrument. (Recall, 
however, that the minimum velocity and temperature 
uncertainties do not occur for the same values of no and 

«f.). 

Experiment 

The optical system described above and shown 
in figure 2 (with the laser output power equal 5 W) was 
used to measure the flow in the potential core of a 
subsonic free jet (15.9 mm exit diameter). The 
laboratory compressed air supply used for the jet was 
passed through a 0.2 pm filter to remove particles. As 
shown in figure 2, the jet axis bisected the incident and 
scattered light directions. With this configuration, the 
measured velocity component was the jet axial 
component. Since the flow velocity was opposite the K 
vector, the frequency shift was negative. Thus the 
Fabry-Perot interferometer fringe diameters of the 
Rayleigh scattered light decreased as the flow velocity 
increased. That is, as the velocity increased, the signal 
from the inner region increased, while the signal from 
the outer region decreased. An electronic pressure 
gauge was used to measure the total pressure in the 
nozzle plenum. These pressure measurements were 
used to calculate the flow velocity and temperature 
using the isentropic flow relations. 

Signals from the three PMT’s were digitized at 
a 50 kHz rate with a digital oscilloscope. For this work, 
the PMT signals were fed directly into the 1 MQ 
oscilloscope inputs. The data record length was 5 sec, 
which corresponded to 250,000 samples on each of the 
three channels. These data were stored on the 
oscilloscope’s hard disk for later processing. 

Data processing 

The data processing scheme was based on the 
model described above. To account for gas density 
changes, the signals from PMT2 and PMT3 were first 
normalized by the signal from PMTl. The next step 
was to find the gain factor that related the measured 
signal voltages to the number of photoelectrons counts 
obtained from the model. To do this, an initial data set 
was taken at a known flow condition, usually zero, or 
close to zero, velocity and a known temperature. The 
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gain factors were obtained by matching the 
photoelectron counts for the two channels from the 
model with the measured values. Note that this 
calibration only requires a measurement at a single flow 
condition. The data taken at other flow conditions were 
then analyzed using the theoretical model to obtain both 
velocity and temperature. Figure 6 shows the grid of 
predicted velocity and temperature as a function of the 
signals from inner and outer regions. The measured 
values of the inner and outer signals (example shown as 
squares on figure 6) were converted to velocity and 
temperature by interpolating between the constant 
temperature and the constant velocity contours. Note 
that the measured values generally follow the model. 
As discussed above, the system was optimized for 
velocity measurement near zero velocity. 

Experimental results 

A typical time record of the digitized signals 
from the three PMT’s is shown in figure 7. The 
bandwidth of the signals was set by the RC time 
constant of the PMT's and associated coaxial cables’ 
capacitance and the 1 M12 oscilloscope input 

impedance (1.12 ms for PMT 2 and PMT 3, determined 
using a pulsed LED as the light source). In figure 8, we 
show an example of a time record obtained while the jet 
total pressure was varied. The inner and outer signals 
are shown along with the velocity obtained from them. 
This illustrates the real-time capability of the technique. 
In figure 9, we show two sets of velocity data taken at 
different fixed jet velocities. The velocity obtained 
from the Rayleigh scattering measurements is shown 
along with the velocity obtained from the measured 
total pressure and isentropic flow relations. The flow 
velocity range was from zero to about -100 m/sec 
(negative because the frequency shift is negative). The 
point used to obtain the calibration is shown on each 
set. It should be noted that the total pressure 
measurement only had a resolution of 0.01 psia, thus the 
isentropic velocities could have a larger error at low 
velocities. These results are based on the average time 
record over the 5 second time interval. The first set 
shows good agreement, including the last point, which 
was taken at the same flow condition as the first point. 
(The first point was used as the calibration point.) The 
second set of measurements shows poor agreement if 
the first point is used as the calibration point. If the last 
point, however, is used as the calibration point, better 
agreement is obtained. The disagreement at the low 
velocity may be a result of the inaccuracies in the 
velocity obtained from the total pressure at low 
velocity. 


Concluding remarks 

Molecular Rayleigh scattering appears to be a 
feasible technique for the dynamic measurement of gas 
density, velocity, and temperature. There are, however, 
several requirements for obtaining successful 
measurements. First, a high power laser source is 
needed to generate a sufficiently large Rayleigh 
scattering signal to overcome the inherent statistical 
photon noise. Second, the maximum amount of the 
Rayleigh scattered light must be collected and detected, 
which means fast collection optics coupled with large 
diameter Fabry-Perot mirrors should be used with low- 
noise, high quantum efficiency photodetectors. The 
Fabry-Perot interferometer is extremely sensitive to any 
outside disturbances, including acoustic noise, 
mechanical vibration, and temperature changes. It is 
highly desirable to locate the Fabry-Perot in a quiet, 
stable environment, with the Rayleigh scattered light 
transmitted through fiber optics from the test facility. 
An automatic system for maintaining the interferometer 
alignment is also essential. 

Two requirements must be met for the 
successful implementation of the technique described. 
One is the elimination of particles in the flow. The 
other is preventing any stray scattered laser light from 
reaching the photodetectors. These requirements will 
probably restrict this technique to unconfined flows 
unless extraordinary measures are taken. It should be 
noted that particle and stray light can be handled if 
more details of the spectrum are measured, which 
means the use of additional detectors. The system 
described here was designed to use the minimum 
number of detectors. 
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Fig. 1 - Comparison of Gaussian and Tenti S6 Rayleigh scattering spectra; y = 0.79. 



Fig. 2 - Setup for Dynamic Rayleigh scattering experiment. 
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Fig. 3 - Fabry-Perot interferometer fringe pattern with mask showing two active areas; n„ — 0.087, «/ 0.039, 

0.341; left image shows full fringe for extended source; right image shows fringe for line source used in experiment. 
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Fig. 4 - Lower bounds for uncertainties for velocity and temperature as a function of rio and nf. (a) and (b) show 
contours of velocity uncertainty in m/sec; (c) and (d) show contours of temperature uncertainty in K; X on (a) and 
(c) at values of no and nf. used for experiment 
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Fig. 5 ~ Uncertainty for velocity (a) and temperature (b) as a function of velocity and temperature for = 0.087 and 
nf= 0.039. 
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INNER COUNTS 

Fig. 6 - Plot of signal from inner and outer active regions as a function of velocity and temperature for n, 
and «/= 0.039. Squares mark location of experimental measurements. 



Sannples 

Fig. 7 - Typical PMT signals (density is divided by 10), sampling rate = 50 kHz. 
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Fig. 8 - Time history of signal from two PMT’s and velocity obtained from these signals. Velocity was manually 
changed over the time period. 
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Fig. 9 - Comparison between time average velocity obtained from isentropic flow relations and from Rayleigh 
scattering measurements. 
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